Voltage induced control and magnetoresistance of noncollinear frustrated magnets 
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Noncollinear frustrated magnets are proposed as a new class of spintronic materials with high 
magnetoresistance which can be controlled with relatively small applied voltages. It is demonstrated 
that their magnetic configuration strongly depends on position of the Fermi energy and applied 
voltage. The voltage induced control of noncollinear frustrated materials (VCFM) can be seen as a 
way to intrinsic control of colossal magnetoresistance (CMR) and is the bulk material counterpart of 
spin transfer torque concept used to control giant magnetoresistance in layered spin-valve structures. 

PACS numbers: 72.15.Eb, 72.15.Gd, 75.10.Lp, 75.10.Jm, 75.47.-m, 75.47.Pq 



The discovery of giant magnetoresistance (GMR) 0,0 
in magnetic multilayered structures has generated a new 
field of spin-based electronics [1,0], or spintronics, which 
combines two traditional fields of physics: magnetism 
and electronics. A spin- valve concept [5J used in GMR 
structures allows controlling the magnetic configuration 
of its ferromagnetic layers by (i) application of relatively 
small magnetic fields or (ii) passing spin polarized cur- 
rents using spin transfer torque (STT) [1, 0j- These fac- 
tors made them ideal systems for spintronic applications 
such as magnetic random access memories (MRAM) and 
magnetic field sensors used in read heads. 

The advent of GMR has considerably increased an in- 
terest in related phenomenon in bulk materials, colossal 
magnetoresistance (CMR) [8j, |9j, which is several orders 
higher than GMR and unlike the latter, can be viewed 
as an "intrinsic" property of material itself. To date, 
the CMR is typically observed in certain manganite com- 
pounds with the bulk magnetic configuration controlled 
by applying the magnetic field (similar to method (i) 
mentioned above for spin valves) but requires character- 
istic magnetic fields of several Tesla [§] . Such high fields 
make them inappropriate for use in spintronic applica- 
tions where appropriate scale should be about Oersteds. 
However, one may expect the possibility of controlling 
the intrinsic magnetic configuration of the bulk materi- 
als (and thus of CMR) by passing spin polarized currents 
through them similar to STT mechanism (ii) mentioned 
above for spin- valves. Since the STT in the latter orig- 
inates from non-collinearity of their adjacent magneti- 
zations, the same requirement should hold for the bulk 
materials. 

Here we promote for the first time magnetically frus- 
trated bulk materials as a new paradigm for spintronic 
applications with high magnetoresistance which can be 
controlled with relatively small applied voltages and does 
not require spin polarized currents. This novel phe- 
nomenon may be viewed as the "bulk" counterpart of 
STT in layered spintronic structures (spin valves) and 
may represent a crucial interest both from applications 
and fundamental points of view. Below we demonstrate 
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FIG. 1: (a) "4in" and (b) "2in2out" configurations for square 
lattice 



that the magnetic configuration of the bulk frustrated 
material is changed under applied voltage leading to the 
strong variation of its conductance. 

The key mechanism at stake is a somewhat microscopic 
equivalent of spin torque, allowing for local spin flips 
within the microscopic spin texture. Those local moves 
requires non collinearity of the spin configuration, which 
would otherwise be insensitive to the current. Further- 
more, in order to avoid a continuous response of the spin 
configuration, we require the ground states to be well 
defined minima, such that the switching from one to the 
other looks like a logical operation. These two aspects are 
well taken into account as soon as the underlying local- 
ized moments possess strong multiaxial anisotropics, in 
order to behave bitwise like and fulfill the non collinear- 
ity criterion. Several Rare-Earth-transition metal inter- 
metallics systems are relevant candidates : in these sys- 
tems frustration arises due to the competition of crystal 
field anisotropy, exchange and quadrupolar interactions 
(for a review see Ref. [10(). These intermetallics systems 
often show a non-collinear magnetic structure: this is 
the case for example of TbGa9 jlll] ■ HoGeafll] or Ura- 
nium compounds [13j|. Naturally, pyrochlores in which 
frustration is due to the crystal structure are also rele- 
vant candidates, as it is well known that in such exotic 



2 



v -7 



CD 



-8 



>"6 
CD 



b} 



-7 



1 1 


1 ' 1 ' 1 ' 1 

Ain / 
*+i 1 1 / 

2in2out /, 




\/ " E=-2tW 


H 1— 

■(b) 


— 1 1 1 1 1 1 h 

Ain / 

*+ii i / 

2in2out / 

\ / x 


1 


s 'e=-2.5 eV 

F 

i i . i . i 



0.0 0.5 1.0 1.5 2.0 

Voltage (V) 



FIG. 2: (a) and (b): total energy as a function of applied 
voltage for Ef~-2 eV and Ep=-2.5 eV, respectively, shown 
by arrows in Fig. [3] 



systems, non collinear low temperature magnetic phases 
are often stabilized [14[, among those, few metallic com- 
pounds have been identified like Pr 2 Ir 2 7 [H[Ii 

and the 

family of the pyrochlore Molybdates, R2M02O7 ljj- 



For argument's sake, we have used a two-dimensional 
square lattice model in order to focus on the driving 
mechanism that is at stake and to describe its qual- 
itative behavior. Realistic realizations would involve 
more complicated structures, like the one in the 3-D py- 
rochlores 14J. We consider a square lattice of classical 
localized moments Sj with strong local on-site uniaxial 
anisotropy (Do) along square diagonals (represented by 
unit vector n,) and coupled through intersite exchange 



lij. In addition, moments Sj are coupled with conduction 
electrons through the local exchange (Jo) described here 
quantum mechanically using tight-binding model. Unlike 
collinear systems, non-collinearity provided by magnetic 
frustration is a key ingredient for switching phenomenon 
proposed here since its origin is due to STT mechanism 
acting locally by conduction electrons on localized mo- 



ments Sj. The Hamiltonian of the system has the form: 

* E + h.c) -J J2 t ■ Si) ef , (i) 

where and Jo are the intersite and the local exchange 
constants respectively, Dq is the uniaxial anisotropy con- 
stant, t is the hopping integral between two neighboring 
sites, c| CT and cf are the creation and annihilation oper- 
ators of the conduction electron with the spin a on site 
% and a a p is the vector of Pauli matrices. For the cho- 
sen model, nearest neighbor interactions I± are irrelevant 
(because nearest neighbor spins are always orthogonal) so 
that only the second nearest neighbor one I2 is taken into 
account. We emphasize that dealing with longer ranges 
interaction does not affect further results reported here 
since only energy differences between magnetic configu- 
rations matter. 

In the limit of low temperature and Dq — > 00 the lo- 
calized moments Si are strictly collinear with rij (Sj||rij) 
and the total energy of the system within a constant shift 
becomes 

Etot = ~h S * ■ S 3 + Tr \Hp\ = 

= -hJ2s i -s j -±-j2f EG U E ) dE ^ ( 2 ) 

i,j 3 

where G < and p represent the "lesser" non-equilibrium 
Green function and density matrix respectively [18l | , and 
the Green function indices include spin index as well. 

We adopt a conventional transport approach where the 
system is considered to consist of three regions: semi- 
infinite left (L) and right (R) regions connected to the 
middle scattering region (M) which is assumed to be 
very long. Despite of all three regions being identical, 
such a subdivision is chosen in order to consider all non- 
equilibrium processes occuring in the scattering region 
M while the left and right regions are in thermodynam- 
ical equilibrium described with Fermi-Dirac distribution 
functions f L ( R ) = .f{E-p L{K) ), where p L ( R ) is the chem- 
ical potential in the left (right) region. The "lesser" Green 
function can be written as 



(3) 



where ip^ is the single-electron wave function on site 
j £ M incident from the left (right). It is straightforward 
from the definition of retarded and advanced Green func- 
tions that 
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(4) 



Since the middle region is assumed to be very long, the fi- 
nite potential drop eV — fiL—^R results into infinitesmal 
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change from one unit cell to another within the middle 
region M. In this case the reflections from the L\M and 
R\M boundaries are negligible and one can write for the 
left(right) wave functions V'.f = c L ^ exp (+(— )k x aj), 
where k x is the wave vector along x-axis and a is the lat- 
tice constant. One can show that the expression for the 
charge current which is proportional to (/x — fit) uM * s 
satisfied when \ipf\ 2 — \ipf\ 2 - Indeed, the charge current 
can be expressed using the lesser Green function as : 



Using Eq. ([3]) and aforementioned wave functions it yields 
the following form: 

I~f**k.a[f L \c L \>-f Il \c ll \>]dE 

which is proportional to — /fl) only when \cl\ 2 = \cr\ 2 
(and consequently \ipf\ 2 — \ipf\ 2 ). 

Then, it straightforwardly follows from Eq. |3]) and 
Eq. (4} that 



GfJE)=i(f L + f R )Z(GUE)) 



(5) 



where 5 notation for imaginary part is used. The total 
energy given by expression ^ of the system per a unit 
cell yields: 



Etot — —I2 J]] Si ■ Sj — 

^J(fL + In) E (E))EdE, (6) 

3 

where G^j(E) is calculated as jth diagonal element of 
the matrix [E — H + iS]^ 1 . Thus, the magnetic config- 
uration corresponding to the minimum of the total en- 
ergy is resulting from the interplay of intersite exchange 
interaction between localized moments on the next near- 
est neighbour sites (the first term in Eq. ©) and the 
exchange interaction between the localized moment and 
the conduction electron described by the second term in 
Eq. ©. 

For simplicity, among the possible magnetic configura- 
tions of the system (Sj||nj) we select "4in" and "2in2out" 
ones which are the two configurations with the largest en- 
ergy separation in this example (see Fig. |T|) . We choose 



21] 



for parameters Jo = 2 eV, t = 1 eV and I2 = O.lt 

In Fig. [5] we show the dependence of the total energy 
for the two aforementioned selected configurations as a 
function of applied voltage. First, one can note that E to t 
strongly depends on the Fermi level position which is 
used to represent different materials. Indeed, when the 
system is in equilibrium, its "2in2out" ("4in") state is en- 
ergetically more favorable for Ep = —2 eV (—2.5 eV) as 
shown in Fig. [21(a) and (b), respectively. Furthemore, the 
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FIG. 3: Total density of states for "4in" (black solid) and 
" 2in2out" (red dashed) configurations represented in Fig. [T] 



strong variation of the total energy as a function of ap- 
plied voltage causes the system to switch from an equilib- 
rium state to another one at certain critical voltage (see 
Fig. [2]). Moreover, further increase of the applied volt- 
age may again reverse the state of the system (Fig. [^t). 
This demonstrates that the magnetic configuration of the 
system can be controlled by the applied voltage. 

The mechanism of these dependences may be under- 
stood from the corresponding total density of states 
(DOS) (Hi for "2in2out" and "4in" configurations repre- 
sented in Fig.[3]where for convenience purposes only neg- 
ative energy range is shown since DOS(-E)=DOS(— E). 
The DOS for both magnetic states have sharp peaks and 
band gaps causing strong dependence of the preferable 
configuration on the Fermi level Ep both in and out of 
equilibrium. Indeed, as on can see from the second term 
111 Eq. ©, the total energy is defined by the sum of two 
products of the DOS with E integrated till hl and n R . 
In the absence of applied voltage, E tot is calculated with 
Mi = l l R = Efi i-e. defined by the position of the Fermi 
level indicated by black arrows in Fig. [3] corresponding 
to two values —2 eV and —2.5 eV used in Fig. [2a) and 
(b), respectively. 

One might expect that the conductances of two afore- 
mentioned magnetic configurations may strongly dif- 
fer due to significant differences in their DOS, thereby 
leading to high magnetoresistance values defined as 
(o^w _ w here &jg and a 2 ™ 2out are the 

conductances for "4in" and "2in2out" magnetic config- 
urations, respectively. We calculate the linear response 
conductance using the Kubo formula [22j 



-Tr 



v x S(E - H)v x S(E - H) 



where v x is the electron velocity operator along the ap- 
plied voltage direction, L — > 00 is the length of the system 
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tremely high values and their magnetic configuration 
may be controlled by applied voltage. The proposed 
phenomenon is the bulk material analog of spin trans- 
fer torque used in layered spin-valve structures. The 
model chosen here used to prove the concept, in combina- 
tion with inverse approach of designing the crystal struc- 
ture from predefined electronic structure properties 24j . 
opens the way towards experimental realization of the 
proposed phenomenon. 

This work was supported by Chair of excellence Pro- 
gram of the Nanosciences Foundation in Grenoble. 



FIG. 4: Conductance as a function of the Fermi level. 



and 

S(E — H) = --%(G r (E)) = — (G r (E) ~ G a (E)) 

TT 7T 

The calculated conductance of the system for "4in" 
and " 2in2out" magnetic configurations are represented in 
Fig. [U As one can see, both a^™ anc i a 2in2out strongly 
depend on position of the Fermi energy Ep and strongly 
correlate with DOS picture (see Fig. [3|. For the two 
values of the Fermi energy used above, the MR ratio is 
found to be about 336 % for the case when Ep — —2.5 eV 
and becomes infinite for Ep = —2 eV. 

One can note that similar behavior is expected to oc- 
cur in systems with collinear antiferromagnetic structure: 
depending on the band structure and the band filling, 
the ferromagnetic state can be close in energy with the 
AF one, and an applied voltage could stabilize the fer- 
romagnetic state. However in such a case, there is no 
spin torque since the structure is colinear, and no mech- 
anism for spin reversal (except through spin waves). On 
the contrary, the mechanism proposed here neither re- 
quire any defect, spin wave nor a fine tuning of any kind. 
The noncollinear nature of the magnetic configurations 
considered provides a robust mechanism for microscopic 
magnetic switching. Furthermore, this mechanism does 
not require spin polarized current injection. 

In conclusion, we have shown that magnetoresistance 
of noncollinear frustrated bulk magnets may reach ex- 



[1] M. Baibich et al, Phys. Rev. Lett. 61, 2472 (1988) 

[2] G. Binasch, P. Griinberg, F. Saurenbach and W. Zinn, 

Phys. Rev. B 39, 4828 (1989) 
[3] A. Fert et al, Mat. Sci. Eng. B, 84, 1 (2001) 
[4] S. Wolf, Science, 294, 1488 (2001) 
[5] B. Dieny et al, Phys. Rev. B 43, 1297 (1991) 
[6] J. Slonczewski, JMMM 159, LI (1996) 
[7] L. Berger, Phys. Rev. B 54, 9353 (1996) 
[8] R. von Helmolt, J. Wecker, B. Holzapfel, L. Schultz, and 

K. Samwer, Phys. Rev. Lett. 71, 2331 (1993) 
[9] S. Jin et al, Science, 264, 413 (1994) 
[10] D. Gignoux and D. Schmitt, J.Magn. Magn. Mater, 100, 

99 (1991) and J. Alloys and Comp. 326, 143 (2001). 
[11] I. Auneau et al, Physica B 212, 351 (1995) 
[12] P. Schobinger-Papamantellos et al, J. Phys. : condens 

Matter 20, 195202 (2008) 
[13] P. Burlet et al, Sol. St. Comm. 39, 745 (1981) 
[14] J. Gardner, M. J. P. Gingras, J. E. Greedan, 

arXiv:0906.3661 ; to appear in Rev. Mod. Phys 
[15] S. Nakatsuji et al, Phys. Rev. Lett. 96, 087204 (2006) 
[16] Y. Machida et al, Phys. Rev. Lett 98, 057203 (2007) 
[17] I. Kezsmarki et al., Phys. Rev. Lett. 93, 266401 (2004)) 
[18] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics, 
Course of Theoretical Physics (Pergamon, Oxford, 1981), 
Vol. 10 

[19] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992) 

[20] C. Caroli and R. Combescot and P. Nozieres and D. 

Saint-James, J. Phys. C: Solid St. Phys., 4, 916 (1971) 
[21] I. V.Solovyev, Phys. Rev. B 67, 174406 (2003) 
[22] R. Kubo, J. Phys. Soc. Jpn, 12, 570 (1957) 
[23] The Density of States is equal to -I/tt^j QGfj(E), 

where summation includes all atoms in the unit cell and 

spins. 

[24] A. Franceschetti and A. Zunger, Nature 402, 60 (1999) 



